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The pseudo-rapidity dependence of the elliptic flow at various excitation energies mea- 
sured by the PHOBOS Collaboration in Au+Au collisions at RHIC is one of the surprising 
results that has not been explained before in terms of hydrodynamical models. Here we 
show that these data are in agreement with theoretical predictions based on perfect fluid 
hydrodynamics. We also show that these PHOBOS data satisfy the universal scaling 
relation predicted by the Buda-Lund hydrodynamical model, based on exact solutions of 
perfect fluid hydrodynamics. 

1. Introduction 

One of the unexpected results from experiments at the Relativistic Heavy Ion Collider 
(RHIC) is the relatively strong second harmonic moment of the transverse momentum 
distribution, referred to as the elliptic flow. Measurements of the elliptic flow by the 
PHENIX, PHOBOS and STAR collaborations (see refs. [E 121 13 EI) reveal rich details in 
terms of its dependence on particle type, transverse and longitudinal momentum variables, 
on the centrality and the bombarding energy of the collision. In the soft transverse 
momentum region, these measurements at mid-rapidity are reasonably well described 
by hydrodynamical models [EJIH]- However, the dependence of the elliptic flow on the 
longitudinal momentum variable pseudo-rapidity and its excitation function has resisted 
descriptions in terms of hydrodynamical models. 

Here we show that these data are consistent with the theoretical and analytic predictions 
that are based on perfect fluid hydrodynamics. 



2. Results 

Our tool in describing the pseudorapidity-dependent elliptic flow is the Buda-Lund 
hydrodynamical model. The Buda-Lund hydro model [Cj is successful in describing the 
BRAHMS, PHENIX, PHOBOS and STAR data on identified single particle spectra and 
the transverse mass dependent Bose-Einstein or HBT radii as well as the pseudorapidity 
distribution of charged particles in Au + Au collisions both at ^/snn = 130 GeV [ |H] 
and at a/snn = 200 GeV [ E] . However the elliptic flow would be zero in an axially 
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symmetric case, so we developed the ellipsoidal generalization of the model that describes 
an expanding ellipsoid with principal axes X, Y and Z. Their derivatives with respect to 
proper-time (expansion rates) are denoted by X, Y and Z. 

The generalization goes back to the original one, if the transverse directed principal 
axes of the ellipsoid are equal, ie X = Y (and also X — Y). 

One can define a deviation from axial symmetric flow, a momentum-space eccentricity: 



X 2 — Y 2 

e p = -. ^. 1 

P X 2 + Y 2 

The exact analytic solutions of hydrodynamics (see ref. [QUI El 02] )> which form the 
basis of the Buda-Lund hydro model, develop Hubble-flow for late times, ie X ^ T ~><x> Xt, 
so the momentum-space eccentricity e p nearly equals space-time eccentricity e. Hence, in 
this paper we extract space-time eccentricity (e) and average transverse flow (ut) from 
the data, instead of X and Y. 

In the time dependent hydrodynamical solutions, these values evolve in time, however, 
it was show in ref. [EI] that X and Y, and so e and u t become constants of the motion 
in the late stages of the expansion. 

The result for the elliptic flow (under certain conditions detailed in ref [ I14j) is the 
following simple universal scaling law: 

The model predicts an universal scaling: every v-i measurement is predicted to fall on 
the same universal scaling curve h/Io when plotted against w. 

This means, that t>2 depends on any physical parameter (transverse or longitudinal 
momentum, center of mass energy, centrality, type of the colliding nucleus etc.) only 
trough the scaling paremeter w. 

Here w is the scaling variable, defined by 

W ^(±-±-\, (3, 



and 



T*,x = T + m t X 2 °_ 2 , (4) 
T + m t a z 



= T + m t Y 2 ( r,) 

and 



m t = m t coshes - y). (6) 

Here a = (AT/T) t represents the temperature gradient in the transverse direction, at 
the freeze-out, m t is the transverse mass, T the central temperature at the freeze-out, 
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while 7] s is the space-time rapidity of the saddle-point (point of maximal emittivity). This 
saddlepoint depends on the rapidity, the longitudinal expansion, the transverse mass and 
on the central freeze-out temperature: 

where y = 0.5 log( ^^ z ) is the rapidity and An represents the elongation of the source 
expressed in units of space-time rapidity. See more details in ref. [Ej. 

Eq. |21 depends, for a given centrality class, on rapidity y and transverse mass m t . Before 
comparing our result to the ^2(77) data of PHOBOS, we thus performed a saddle point 
integration in the transverse momentum variable to end up with a formula that can be 
directly fitted to v 2 (v) of PHOBOS. 

The fitting package is available at ref. [E3, more about the results (eg. contour plots) 
are available at ref. [I16|. 

We have found that the essential fit parameters are e and An, and the quality of the 
fit is insensitive to the precise value of T , a, u t and R g . These parameters dominate the 
azimuthal-averaged single particle spectra as well as the HBT (Bose-Einstein) radii, how- 
ever they only marginally influence v%. In a broad region their precise value is irrelevant 
and does not influence the confidence level of the t>2(?7) fits. Hence we have fixed their 
values as given in the caption of tabled We also removed points with large rapidity from 
the fits in case of lower center of mass energies. 

Fits to PHOBOS data of ref. [[J are shown on the left panel of fig. ^ The right panel 
of fig. [H demonstrates that these data points follow the theoretically predicted predicted 
scaling law. 





19.6 GeV 


62.4 GeV 


130 GeV 


200 GeV 


e 

An 

X 2 /NDF 
conf. level 


0.294 ± 0.029 
1.70 ± 0.25 
1.8418/11 
99.995% 


0.349 ± 0.008 
2.16 ± 0.05 
20.1388/13 
21.4036% 


0.376 ± 0.005 
2.46 ± 0.04 
34.7935/15 
1.00341% 


0.394 ± 0.006 
2.56 ± 0.04 
27.4865/15 
7.03121% 



Table 1 

Results of fits to PHOBOS data of ref. [GJ. Remaining parameters were fixed as follows: 
T = 175 MeV, a = 1.19 and u t = 1.64. 



3. Conclusions 

In summary, we have shown that the excitation function of the pseudorapidity depen- 
dence of the elliptic flow in Au+Au collisions is well described with the formulas that are 
predicted by the Buda-Lund type of hydrodynamical calculations. 

We have provided a quantitative proof of the validity of the perfect fluid picture of soft 
particle production in Au+Au collisions at RHIC but also show here that this perfect 
fluid extends far away from mid-rapidity. 
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Buda-Lund fits to PHOBOS v (n) data 



PHOBOS v. versus scaling variable w 
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Figure 1. On the left: PHOBOS data on the pseudorapidity dependence of the elliptic 
flow [ 1 , at various center of mass energies, with Buda-Lund fits. On the right: elliptic 
flow versus scaling variable w is plotted. The data points show the predicted scaling 
behavior 



The universal scaling of PHOBOS v 2 (r]), expressed by eq. 121 and illustrated by fig. [T] 
provides a successful quantitative as well as qualitative test for the appearence of a perfect 
fluid in Au+Au collisions at various colliding energies at RHIC. 
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